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Abstract 

A search for the pair production of a heavy, narrow resonance decaying into two 
jets has been performed using events collected in a/s = 7 TeV pp collisions with the 
CMS detector at the LHC. The data sample corresponds to an integrated luminosity 
of 5.0fb"\ Events are selected with at least four jets and two dijet combinations with 
similar dijet mass. No resonances are found in the dijet mass spectrum. The upper 
limit at 95% confidence level on the product of the resonance pair production cross 
section, the branching fractions into dijets, and the acceptance varies from 0.22 to 
0.005 pb, for resonance masses between 250 and 1200 GeV. Pair-produced colorons 
decaying into qq are excluded for coloron masses between 250 and 740 GeV. 
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The high center-of-mass energy provided by the Large Hadron Collider (LHC) offers oppor- 
tunities to search for physics beyond the standard model (SM) and in particular to search for 
new strongly interacting particles. Searches for new resonances in the dijet mass spectrum for 
the jets with the highest transverse momentum (pj) have been performed at both the Teva- 
tron [IJ and the LHC ||2HZI- These searches were not optimized for pair-produced particles, 
which are predicted by some models ||8l- [T0l . The ATLAS experiment has performed a search 
for pair production of dijet resonances in four-jet events fTll.'T2'| and excludes a scalar gluon at 
95% confidence level (CL) in the mass range between 100 and 287 GeV. 

In this Letter, we report the results of a search for the pair production of a narrow resonance, 
which decays into two jets, using the dijet mass spectrum in four-jet final states, measured with 
the Compact Muon Solenoid (CMS) detector in pp collisions at the LHC at i/s = 7 TeV. This 
search focuses on the high mass range between 250 and 1200 GeV. We define the common dijet 
mass as the average of the invariant masses of the two dijets with the smallest difference in 
invariant mass, rejecting events where the difference exceeds 15% of the average value. Mod- 
els IHI predicting pair production through gg interactions of color-octet vectors, also called 
colorons (C), and color-octet scalars (Ss) would give rise to such a signature. In addition to qq 
decays, the coloron can also decay into a pair of color-octet scalars if the mass of the coloron 
is greater than twice the mass of the two scalars. Pair-produced colorons decaying to quark- 
antiquark pairs (gg — > CC — > qqqq) fS] are used as the benchmark signal in this analysis, 
although we separately consider the possibility of decays to SsSg. As a third possibility, we 
consider an R-parity violating SUSY model 1(131 [r4 | in which pair-produced top squarks (stops) 
each decay to qq, as this leads to a similar final state. 

The CMS detector is a multipurpose apparatus and is described in detail in Ref. IITSll . The CMS 
coordinate system has its origin at the center of the detector, the z axis along the direction of 
the counterclockwise circulating proton beam, the y axis normal to the LHC plane pointing 
vertically upward, and the x axis radially inward toward the center of the LHC ring. We de- 
fine (p to be the azimuthal angle, 6 the polar angle, and f] = — ln[tan(0/2)] the pseudorapidity. 
The central feature of the CMS apparatus is a superconducting solenoid with a 6 m internal 
diameter, operating at a central field strength of 3.8 T. Within the field volume are, in order of 
increasing radius, a silicon pixel and strip tracker, a high-granularity PbW04 crystal electro- 
magnetic calorimeter (ECAL), and a brass /scintillator hadron calorimeter (HCAL). All three 
systems have both barrel and endcap components, with \rj\ = 2.5 (3.0) the outermost extent 
of the tracker (calorimeters). The ECAL and HCAL cells are grouped into towers, projecting 
radially outward from the origin. Outside of the field volume an iron/ quartz-fiber hadron cal- 
orimeter covers the forward region (3 < |//| < 5). A muon system encloses the central and 
endcap calorimeters out to | 1 = 2.4. 

The data sample used for this analysis was collected in 2011 and corresponds to an integrated 
luminosity of 5.0 fb^^. The triggers used for the analysis require the presence of at least four 
jets, based on information from the calorimeters. Each jet must have |//| less than 3.0 and pj 
greater than 70 or 80 GeV, depending on the running period. This trigger is 99.5% efficient 
for events with four leading (highest pj) jets, each with a transverse momentum exceeding 
110 GeV. 

For offline reconstruction, we employ the CMS particle-flow algorithm IIT6ll in the region < 
2.5 to reconstruct objects used in jet determinations. This algorithm uses calorimeter infor- 
mation and reconstructed tracks to identify electrons, muons, photons, and both neutral and 
charged hadrons. Jets are reconstructed from particle-flow objects using the anti-Zcj algorithm 
with a distance parameter 0.5 ||T7l[l8| . Jet energy is corrected to account for the nonlinearities 
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and nonuniformities in the response of the calorimeters, as determined from Monte Carlo (MC) 
simulation, test beam, and collision data [|i9i . Additional corrections accounting for the effect 
of multiple pp collisions per bunch crossing are also applied ll20ll2T| . 

We require events to have at least one good primary vertex with a z position within 24 cm 
of the center of the detector and with a transverse distance from the beam spot of less than 
2 cm. A set of jet quality criteria are applied to remove possible instrumental and non-collision 
backgrounds [22 J. All data as well as all simulated signal events passing these selection criteria 
also satisfy standard jet identification requirements [23|. We require that events have at least 
four jets, each with pj > 110 GeV and \r]\ < 2.5. We require the two jets in each possible pair 
to have a separation ARjj = \J (Af/)^ + (Ai^)^ > 0.7. This ensures a negligible overlap between 
the jets. We calculate the dijet-mass combinations from the four leading jets and choose the one 
with the smallest Am/wjavg/ where Am is the mass difference between the two dijets and wJavg 
is their average mass. We require Am < O.lSmavg/ which is approximately three times the dijet 
mass resolution of 4.5%. 

The benchmark signal events are simulated using the MadGraph v5 [24] event generator with 
the CTEQ6L1 parton distribution functions (PDF) |[25l, and PYTHIA v6.4.26 lEEl parton show- 
ering and hadronization. The generated events are further processed through a Geant4 [271 
simulation of the CMS detector. The assumed width of the simulated coloron resonance is 
negligible compared with the experimental resolution. The dominant background arises from 
QCD processes resulting in four or more jets. Studies of this background are performed using 
a sample of simulated QCD events also generated using MadGraph. 

For each dijet we define a quantity A as the difference between the scalar sum of the transverse 
momenta of the two jets in the dijet and the average pair mass in the event: A = Yai=\;A.Vi)i ~ 
mavg- Figure [l] shows the distribution of A versus ffJavg for simulated QCD background events 
as well as for coloron signal events. Due to the selection requirements, we observe a broad 
structure at around ffJavg = 300 GeV from QCD events [28J. To remove this structure, thus 
leaving a smoothly falling dijet mass spectrum, we require A > 25 GeV for each of the two 
dijets in the event. This requirement reduces the QCD background by more than an order of 
magnitude while retaining approximately 25% of the signal. 

Figure |2] shows the paired dijet mass spectrum in data with all the selection criteria applied. 
The observed mass spectrum extends up to 1200 GeV. We obtain a prediction for the QCD 
background by fitting the data to a smooth parametrization: 

dc7 ^ Po(l-mavg/v/s)^l 
dmavg (mavg/v/s)^2+P3ln{m,,g/ys)' 

where ¥q, ¥\, P2, and P3 are free parameters. This functional form has been used in previous 
searches for dijet resonances IH. The fit to the data and the normalized QCD simulation are 
given in Fig.Sby solid and dashed-dotted curves, respectively. The fit has a x.^/'^-O-f of 0.94 
over the full mavg rnass range. Although there is an apparent bias toward positive pull values 
in the low mass region, such a bias would result in the quoted limits being conservative in this 
region. 

The signal shapes from colorons and stops have negligible difference and we use a single pa- 
rameterization for both. It is modeled by the sum of two separate Gaussian functions: one 
Gaussian describes the core and the other the tails, with widths and normalizations determined 
from a fit to simulated signal events at each assumed mass value. The dijet mass resolution de- 
scribed by the rms of the core Gaussian is approximately 4.5%, the dijet mass tail described by 
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Figure 1: The distribution of A = E,=i^2(Pt)! — wiavg versus wJavg for QCD multijet background 
alone (main plot) and background plus coloron signal (insets), for simulated events. We remove 
events with an entry below the horizontal line at 25 GeV. There are two entries per event (one 
for each of the two dijet pairs). 



the rms of the other Gaussian is between 150 and 250 GeV, and the fraction of the core Gaus- 
sian varies from 85% at 300 GeV to 45% at 1000 GeV. The signal acceptance, listed in Table [ij 
varies from 0.4% for a coloron with mass 200 GeV to 12.1% for a coloron with mass 1000 GeV. 
The acceptance for the stop signal is larger than that for the coloron signal because the stop 
production model includes qq interactions and has a different final state angular distribution. 

We search for pair production of dijet resonances by fitting the background parameterization in 
Eq. ^ plus a signal hypothesis at an assumed mass to the data. The signal magnitude is a free 
parameter in the fit for each fixed value of signal mass. We restrict the fit to dijet masses above 
220 GeV in order to avoid the threshold due to the jet px requirement. With this requirement the 
dijet mass spectrum is described by both the simulation and the background parameterization. 
The maximum value of the likelihood obtained from a background-only fit is denoted by Lq 
and the likelihood from a background plus signal fit by Lg. The local significance is defined 
as a/— 2 ln(Lo/Ls). The results of this analysis are shown in Fig. [2] We find that the largest 
fluctuation in the pair-produced dijet mass spectrum occurs for a hypothetical resonance mass 
of 1130 GeV and has a local statistical significance of 2.6c. The global significance is reduced to 
l.lcr after taking into account the trials factor [29| within the full mass range of this search. We 
conclude that there is no evidence for pair-produced narrow dijet resonances in the data. 

Upper limits are placed on the product of the production cross section of the pair-produced res- 
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Figure 2: The average paired dijet mass distribution (black points) in data compared with a 
smooth background fit using Eq. ([T]| (blue solid curve) and the result of the fit of the same func- 
tion (dash-dotted red curve) to QCD simulated data (not shown). Also shown are examples of 
simulation of hypothetical coloron resonances decaying 100% to qq (dashed green curves) with 
masses niQ = 400 and 800 GeV. The bin-by-bin fit pulls are shown below. 



onances, the square of the branching fraction to dijets, and the detector acceptance. The dijet 
mass for the limit is required to be above 250 GeV to ensure a full coverage of the low mass tail 
of the resonance between 220 and 250 GeV. To set upper limits we use a modified-frequentist 
method (CLg) p30','3T|. We fit the signal+background hypothesis to the data, allowing both the 
signal strength and the background parameters free to vary. The sources of systematic uncer- 
tainties consist of a luminosity uncertainty of 2.2% ||32| and a signal acceptance uncertainty of 
10%. The latter is determined by the jet energy scale uncertainty (2.2%) and the jet energy reso- 
lution uncertainty (10%) |19|. The variation in expected signal yields due to PDF uncertainties 
is negligible. The uncertainties in the luminosity, the signal acceptance due to jet energy scale 
and resolution, and the parameters of the background function are all treated as nuisance pa- 
rameters and expressed as log-normal distributions with their central values and uncertainties. 
The observed and expected limits are calculated using the CLg method with a one-sided profile 
likelihood test statistic. 

Figure |3] shows the observed and expected 95% CL limits, the la and la uncertainty bands 
around the expected limits, and predictions from the coloron and SUSY models. The observed 
limit on the product of the resonance pair production cross section, the branching fractions 
into dijets, and the acceptance varies from 0.22 to 0.005 pb for resonance masses between 250 
and 1200 GeV. The limits are generally applicable for pair-produced resonances, each decaying 
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to dijets, and they are compared with calculations for the coloron model HI described above. 
At 95% CL we exclude pair production of colorons with mass mc in the range 250 < nic < 
740 GeV, assuming that colorons have flavor-universal couplings and decay only into qq [lOJ. 
Assuming the branching fraction of colorons into qq is reduced due to competition with a 
C — ?> SsSg channel where mgg = 150 GeV and tan0 = 0.3 (the suppression factor of gluon 
coupling to qq compared with the analogous QCD coupling) [10], we exclude pair production 
in the range 250 < nic < 580 GeV. This analysis is not sensitive to the pair-produced Sg, where 
the color-octet scalars decay exclusively to qq. We also compare the results with those of a 
SUSY model for pair-produced stops, where the stops decay exclusively to qq and R-parity 
is violated ||T3l[T4| . The calculation is done at next-to-leading order (NLO) with next-to-NLO 
corrections |l 33ti37| . 



Table 1: The acceptances for the coloron and stop models after applying all selection criteria. 
Most of the variation in the acceptance as a function of resonance mass is due to the jet pj 
requirement. 
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Figure 3: The observed and expected 95% CL limits on the product of the resonance pair pro- 
duction cross section, the square of the branching fraction to dijets, and the detector acceptance, 
given by the solid and dot-dashed black curves, respectively. The shaded regions indicate the 
Icr and 2(7 bands around the expected limits. Predictions of a coloron model and a SUSY model 
are also shown. 

In summary, a search for pair production of a narrow dijet resonance has been performed with 
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the CMS detector using 5.0 fb"^ of v^= 7TeV pp collisions produced at the LHC. The paired 
dijet mass spectrum is found to be a smooth distribution and is in agreement with the predic- 
tions of the standard model. Upper limits are reported on the product of the production cross 
section, the branching fractions into dijets, and the acceptance of a pair-produced dijet reso- 
nance having a width negligible compared with the experimental resolution. At 95% CL, the 
pair production of colorons is excluded for coloron masses between 250 and 740 GeV assuming 
that a coloron decays 100% into qq, or between 250 and 580 GeV assuming that coloron decays 
into qq compete with decays into SsSs- The search significantly extends previous results [12]. 
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